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ABSTRACT 

The  Marsh-Schulkin  (M-S),  or  Colossus,  model  of  acoustic  transmission  loss  in 
shallow  water  is  reviewed  in  light  of  new  information  and  techniques  that  have  emerged 
since  its  introduction  in  1962.  The  M-S  model  is  a  semiempirical  mode)  based  on  exten¬ 
sive  measurements  taken  off  the  Atlantic  Coast.  It  uses  several  concepts:  (1)  refractive 
cycle,  or  skip  distance,  (2)  deflection  of  energy  into  the  bottom  at  high  angles  by  scatter¬ 
ing  from  the  sea  surface,  and  (3)  a  simplified  Rayleigh  two-fluid  model  of  the  bottom  for 
sand  or  mud  sediments.  With  a  few  free  parameters,  including  water  depth,  about 
100,000  measurements  from  100  Hz  to  10  kHz  were  fitted  within  stated  error  bounds. 
The  model’s  chief  criticisms  have  been  that  it  could  not  be  adjusted  for  arbitrary  negative 
sound-speed  gradients  (it  uses  the  same  constant  thermocline  in  all  cases),  and  that  it  uses 
empirical  bottom  loss  values.  The  M-S  model  yields  good  predictions  when  applied 
properly  ;  it  is  not  to  be  used  for  all  environments. 

The  M-S  model  is  compared  with  Rogers*  semiempirical  model  based  on  numerical 
calculations  of  the  normal-mode  solutions  to  the  wave  equation  for  propagation  in  shal¬ 
low  water.  (Rogers’  model  for  transmission  loss  involves  a  simple  range-dependent  fit 
with  three  coefficients,  and  allows  a  convenient  sensitivity  analysis  for  the  addition  or 
changing  of  parameters.)  Differences  in  predicted  propagation  losses  are  discussed. 

The  M-S  model  is  also  extended  for  the  purpose  of  treating  arbitrary  negative  and 
bilinear  sound-speed  gradients.  These  extensions  use  new  general  expressions  for  the 
skip  distance,  the  near-field  anomaly,  and  the  reflection  coefficients.  The  reflection 
coefficients  are  calculated  from  the  Morse-Mackenzic  relations  for  loss  per  bounce  using 
the  values  for  bottom-sediment  properties  reported  by  E.L.  Hamilton.  The  extended 
model  and  Rogers’  model  are  found  to  give  about  tlte  same  predictions  when  the  same 
inputs  are  used. 
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I.  INTRODUCTION 

This  document  is  the  first  of  two  reports  due  under  an  expanded  ONR  contract  con¬ 
cerning  wo;‘k  in  shallow-water  acoustics.  A  great  deal  of  research  has  gone  into  shallow- 
water  acoustics  over  the  years,  especially  as  it  relates  to  the  solution  of  Navy  problems. 
When  the  scientific  and  technical  basis  of  previous  work  has  been  hidden  for  many  years 
in  classified  reports  issued  by  individual  Navy  Centers  and  in  the  minds  of  departed  per¬ 
sonnel  whose  knowledge  has  become  unavailable,  an  information  and  communications 
gap  often  develops.  This  situation  gives  rise  to  needless  duplication  of  effort  and  the  loss 
of  much  important  information  gathered  in  previous  expensive  and  time-consuming 
measurement  programs. 

Such  is  the  case  for  the  Colossus  II  Shallow-Water  Acoustics  Propagation  Studies, 
the  results  of  which  were  detailed  25  years  ago  in  a  classified  document.  Project 
Colossus  n  was  established  in  1954  to  investigate  acoustics  in  shallow  water 
(100  fathoms  or  less).  A  portion  of  that  program  was  devoted  to  a  study  of  underwater 
acoustic  transmission  loss  in  the  frequency  range  of  1(X)  to  3000  Hz.  A  shorter, 
unclassified  report'  and  a  letter  to  the  editor  of  the  Journal  of  the  Acoustical  Society  of 
America^  were  also  prepared. 

The  letter  pi’.blished  in  the  open  literature^  contained  the  useful  transmission-loss 
model  summarizing  almost  ten  years  of  theory,  analysis,  and  at-sea  measurements.  The 
compact  expressions,  including  an  error  table,  were  based  on  about  100,000  observations 
of  propagation  loss  —  also  including  the  results  of  a  far  more  extensive  semiempirical 
data  biise,  Project  AMOS,^  which  contained,  among  other  results,  surface-loss  data  from 
192  acoustic  stations  occupied  mostly  in  the  deep  water  of  the  North  Atlantic  and  the 
Mediterranean  during  the  years  1949  to  1953.  The  frequency  range  covered  by  AMOS 
was  2-25  kHz. 

With  time  our  use  of  the  word  "scmiempiricar*  to  describe  the  resultant  expressions 
gave  rise  to  misconceptions  and  was  reduced  by  others  to  the  word  "empirical."^ 

In  retrospect,  not  allowing  for  arbitrary  negative  gradients  in  tlie  water  may  have 
been  a  major  drawback  of  the  model.  At  that  time,  however,  so  little  was  known  about 
bottom  properties,  and  shipboard  computer  capability  was  so  primitive,  that  a  more 
dC'iuled  predictive  model  was  premature.  The  emphasis  of  the  model  was  on  predicting 
bottom  loss  along  with  the  coupling  of  surface  scattered  energy  to  the  bottom. 

The  original  document  describing  the  M-S  model*  contained  the  following  caveat 
section,  called  "Other  Modes  of  Propagation,"  which  listed  some  of  the  model’s  limita- 
lions; 
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"OTHER  MODES  OF  PROPAGATION 

The  foregoing  description  applies  to  the  propagation  situation  which  occurs  most  of  the 
time  off  the  East  Coast  of  the  United  States.  Internal  channels  were  observed  so  seldom 
that  they  were  not  considered  an  important  situation  in  our  analysis.  Presumably,  the  sur¬ 
face  and  bottom  do  not  affect  this  pixpagadon  mode  which  is  explicable  m  terms  of 
spreading  and  temperature  absorption  losses.  This  mode  is  known  to  be  important  in  cer¬ 
tain  localities  such  as  the  Scotian  Shelf.  Very  little  effect  was  found  with  respect  to 
source  and  receiver  depths.  [Internal  channels  have  a  large  depth  dependence.] 

No  attempt  was  made  to  cany  the  analysis  below  100  cps.  This  frequency  has  a 
wavelengA  about  a  quarter  of  the  water  depth  of  200  feet  typical  of  these  measurements. 
Undoubtedly  a  ray  picture  does  not  apply  here  and  nc^^  mode  theory  is  required  to 
explain  the  observations.  Another  phenomenon  associated  with  low  frequencies  which 
was  not  treated  here  was  propagation  from  one  point  to  another  by  way  of  the  sub-bottom 
or  the  seismic  mode.  The  existence  of  this  path  has  definitely  been  established,  and  can 
be  of  importance  at  longer  ranges  and  lower  fluencies. 

This  analysis  is  not  applicable  to  bottoms  with  sustained  slopes  in  which  a  ray  picture 
accounting  for  progressively  changing  limiting  angles  must  be  employed.  The  shallow- 
water  analysis  applies  up  to  depths  of  100  fathoms,  which  occur  at  about  100  miles  from 
shore. 

In  the  transition  region  between  shallow  water  and  deep  water,  it  is  known  that  propaga¬ 
tion  may  be  described  quite  well  in  terms  of  a  detailed  ray  analjsis." 

When  tested  under  appropriate  situations  tlie  model  has  performed  well.  It  was  not 
meant  to  be  used  in  isovelocity  situations  where  the  Pekeris  normal-mode  model^’^  or  the 
Weston  ray  modef  apply.  Even  under  strong  negative  gradients,  however,  the  M-S 
model  was  attributed  to  give  the  right  answer  for  the  wrong  reason,®  i.e,  because  the 
average  negative  gradient  which  was  used  for  all  cases  still  gave  predictions  within  stated 
error  bounds  (see  Table  V  in  Section  III.D).  In  this  report,  the  M-S  model  has  been 
extended  to  include  arbitrary  negative  and  bilinear  gradients. 

In  Section  II,  we  review  the  M-S  model  and  describe  its  three  main  parameters  — 
skip  distance,  the  near-field  anomaly  correction,  and  the  effective  attenuation  coefficient. 
Section  in  describes  the  Onboard  Prediction  Model  for  Propagation  Loss  developed  by 
P.H.  Rogers  of  the  Naval  Research  Laboratory  (NRL)  and  compares  its  results  with  those 
of  the  M-S  model.  Sections  IV  and  V,  respectively,  describe  how  the  M-S  model  can  be 
extended  to  include  arbitrary  negative  and  bilinear  gradients  and  compare  the  results  of 
the  extended  model  with  those  of  Rogers’  model  under  negative  gradient  circumstances. 
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n.  MARSH-SCHULKIN  SHALLOW-WATER  PROPAGATION-LOSS  MODEL 

Marsh  and  Schulkin'"*  were  among  the  first  to  show  the  importance  of  refractive 
cycles,  or  skip  distances,  for  acoustic  transmission  in  both  deep  and  shallow  water.  In 
shallow  water,  propagation  at  extended  range  and  moderate  frequencies  is  supported  by 
repeated  bottom  and  sixrface  reflections,  almost  regardless  of  the  thermal  conditions. 
Contact  of  rays  with  the  rough  sea  surface  causes  the  scattering  of  energy  at  high  grazing 
angles  and  consequent  loss  into  the  bottom.  Thus,  there  is  a  strong  surface-bottom  cou¬ 
pling  such  that  the  propagation  losses  are  controlled  by  the  number  of  contacts  of  rays 
with  both  surfaces.  The  thermal  structure  of  the  water  affects  propagation  through  its 
influence  on  the  skip  distance  and  the  number  of  surface  and  bottom  contacts.  Determin¬ 
ing  the  skip  distance  and  hence  the  number  of  bottom  bounces  over  an  acoustic  path  is 
important  for  finding  the  losses  suffered  at  the  bottom.  The  vehicle  used  by  Marsh- 
Schulkin  for  establishing  skip  distances  for  both  surface  and  bottom  contacts  was  a 
bilinear  sound-speed  profile  with  a  variable  depth  for  the  positive-gradient  layer. 

Along  with  the  concepts  of  skip  distance  and  bottom  loss,  the  Colossus  model,  also 
known  as  the  Marsh-Schulkin  (M-S)  model,  used  (1)  the  AMOS  results  for  a  deep-water 
isothermal  surface  duct  but  with  the  average  thcrmocline  appropriate  for  shallow  water, 
and  (2)  measurements  of  tlie  actual  propagation  loss  in  shallow  water  off  the  East  Coast 
as  a  function  of  frequency,  separated  by  bottom  type  (sand  or  mud)  and  by  season.  Two 
other  mechanisms  characteristic  of  shallow  water  processes  were  also  included  in  the 
M-S  model:  (1)  a  "near-field  anomaly"  correction  in  the  direct  radiation  zone  that 
included  the  gain  due  to  multiple  bottom  and  surface  bounces,  and  (2)  an  energy  conser¬ 
vation  rule  was  used  to  establish  the  effective  shallow-water  attenuation  coefficient,  a, , 
which  includes  the  additional  loss  due  to  the  coupling  of  energy  from  the  the  wind- 
roughened  sea  surface  to  the  bottom. 

The  bilinear  gradient  used  in  the  M-S  model  is  composed  of  two  constant,  linear 
segments  drawn  toward  the  surface  and  toward  the  bottom  from  the  depth  of  maximum 
sound  speed  (or  temperature).  The  sound  ray  cycles  have  one  upward  radius  of  curvature 
(positive  sound-speed  gradient)  for  surface  bounces  and  one  downward  radius  of  curva¬ 
ture  (negative  gradient)  for  bottom  bounces.  Based  on  the  depth  of  the  surface  layer  and 
the  water  depth,  a  single  effective  skip  distance  is  formulated.  Multiples  of  this  effective 
skip  distance  are  used  to  define  a  zone  of  direct  ray  paths  (20  log/? ,  where  R  =  range),  a 
zone  of  mode  stripping  (1.5  log/?),  and  a  zone  of  single-mode  control  (10 log/?).  The 
mode-stripping  process  was  found  to  be  complete  at  a  range  equal  to  8// ,  where  H  is  one 
effective  skip  distance. 
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A.  Propagation  Loss  Equations 

In  the  M-S  model,  the  propagation  loss  is  thus  represented  in  terms  of  sea  state 
(wave  height),  bottom  type  (or  bottom  loss,  if  known),  water  depth,  frequency,  and  the 
depth  of  the  positive-gradient  layer.  The  skip  distance  is  used  as  a  reference  to  dehne 
regions  where  wave-front  spreading  follows  square,  three-halves,  and  first-power  laws  as 
a  function  of  range.  If  the  range  R  between  source  and  receiver  is  less  than  or  equal  to 
the  skip  distance  H ,  the  propagation  loss  N  is 

N  =2QlogR +aR +60- ki  dB,  (1) 

where  R  is  the  range  in  kiloyards,  a  is  the  absorption  coefficient  in  seawater  in 
decibels/kiloyard,  and  kj^  is  the  near-field  anomaly.  For  intermediate  ranges, 
H  Sf?  S8//, 


=  15  log/? -ha/? 1^^- lJ-F51og// -F60-4  dB-  (2) 

For  long  ranges,  /?  ^  8/f , 

^  -  •< 

/V  =  lOlog/? +a/? ^  ^  +  10  log// +  64.5  -  dB.  (3) 

H 

to  « 

These  equations  provide  for  the  gradual  transition  from  spherical  spreading  in  the  near 
field  to  cylindrical  spreading  in  the  far  field. 

We  now  discuss  the  parameters  that  are  used  in  the  model. 


B.  Skip  Distance 

It  was  stated  earlier  that  the  M-S  model  used  the  concept  of  skip  distance  for  acous¬ 
tic  propagation  in  shallow  water.  The  skip  distance  is  defined  in  such  a  way  that  it 
represents  the  maximum  range  at  which  rays  first  make  contact  with  either  the  sea  surface 
or  the  bottom  and  thus  specifies  the  near-field  region. 

In  defining  skip  distance,  M-S  used  a  single,  negative  sound-speed  gradient  for  the 
layer  below  the  thermocline  and  a  single,  positive  gradient  for  the  layer  above.  Both 
values  were  based  on  the  average  values  observed  during  the  Colossus  measurements: 
-i-O.OlS  s”'  for  the  upper,  isothermal  layer  and  -0.035  s'*  for  the  lower  layer.  If  we  divide 
the  range  by  the  skip  distance,  we  get  the  effective  number  of  bottom  contacts  that  con¬ 
tribute  to  the  propagation  loss. 
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The  original  M-S  skip  distance  // ,  in  kiloyards,  is  given  by 


L+D 


where  D  is  the  depth  of  the  water,  in  feet,  and  L  is  the  depth  of  the  mixed,  isothermal 
layer,  in  feet. 

The  M-S  skip  distance,  in  kilometers,  is 


r 

1“  • 


where  L  and  D  are  in  meters. 


C.  Near-Field  Anomaly 

Results  of  the  Colossus  measurements  showed  that  at  short  ranges  (i.e.,  in  the  near 
field)  propagation  was  consistently  superior  to  that  predicted  for  inverse  square-law 
spreading.  The  departure  from  square-law  spreading  in  the  near  field  can  be  accounted 
for  by  gains  due  to  boundary  reflections.  The  M-S  model  uses  empirical  boundary  loss 
values  determined  from  measurements. 

Figure  1  is  a  ray  diagram  for  the  direct  radiation  zone.  This  diagram  can  be 
expected  to  apply  in  the  region  that  is  essentially  free  from  refraction.  It  can  be  seen  that 
there  is  one  direct  ray  (labeled  1  in  the  figure)  with  no  boundary  contacts,  one  ray  having 
one  more  surface  contact  than  bottom  contact  (labeled  2),  one  ray  having  one  more  bot¬ 
tom  contact  than  surface  contact  (2').  and  thereafter  two  rays  for  each  order  having  an 


Source 


Receiver 
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equal  number  of  bottom  and  siuface  contacts  (e.g.,  3  and  3')-  The  set  of  rays  is  complete. 
Accordingly,  upper  and  lower  limits  on  the  sound  field  £  can  be  calculated  assuming 
incoherent  addition  of  rays: 

Z  =  /  +  2/  (r,  +  r}  +  ....  +  r,"  r^")  (5) 

+  /(ri,+r,  rj2  +  ....  +  r;-lrfc") 


Z  =  /  +  /(2r,  +  r,  +  r|,)(l  +  r,  +  r}  +  ....  +  r,""*  rj,""*) 


Z  =  / 


j  ^  (2r,  Tf, -Fr^Kl 

(1-r,  r^) 


where 


fsfb  ^  1  . 


I  =  intensity  of  direct  ray 
r,  =  surface  reflection  coefficient 
=  bottom  reflection  coefficient. 


The  term  in  brackets,  K^  u,  is  called  the  near- field  anomaly. 


^L,u  =  ^  ••• 


4-  (l-r,%") 


(6) 


where  n  is  the  number  of  bottom  and  surface  contacts  that  contribute  to  the  field  in  the 
near  zone,  defined  as  the  first  skip  distance,  n  is  limited  by  the  critical  angle.  The  multi¬ 
plier  has  an  upper  limit. 


f(u  = 


and  a  lower  limit, 

^1  -F  (2r,ri,+r,+r^)  j  . 

The  gain,  in  decibels,  due  to  the  near-field  anomaly  is 


^L.u  -  lOlogKi  u  . 


(7) 

(8) 

(9) 
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Reflection  coefflcients  r,  and  rj,  were  believed  to  decrease  rather  rapidly  with 
increasing  grazing  angle  of  incidence.  If  this  were  true,  might  be  a  better  representa¬ 
tion  than  ky  in  the  direct  radiation  zone,  or  near  field.  The  near  field  is  defined  as  ranges 
less  than  or  equal  to  the  direct  path  range,  which  is  taken  as  equal  to  the  skip  distance. 

Only  ki  was  used  in  the  published  results,  but  there  is  reason  to  believe  that  k(j  fits 
the  data  better  under  certain  conditions  (see  Table  lA). 

For  reference,  we  list  in  Table  IB  values  for  ki  and  ky  over  sand  and  mud  bottoms 
for  sea  state  2  conditions  and  four  frequencies. 


Table  I  A.  Median  measured  minus  calculated  tosses  (k^^ky  compared).^ 


Frequency, 

Hz 

3 

9 

Range,  kyd 
30 

60 

90 

Near-Field 

Anomaly 

112 

-3 

-4 

-5 

-3 

-3 

h 

446 

-3 

-5 

-4 

-4 

-4 

1120 

-2 

•1 

-2 

-1 

-7 

2820 

-1 

2 

2 

-6 

-9 

112 

0 

0 

2 

0 

1 

ky 

446 

-1 

0 

0 

-3 

-6 

1120 

-2 

0 

5 

5 

-5 

‘'The  data  set  for  ki^  and  ky  were  somewhat  different  by  geography  and  by 
bottom  type.  Thus,  compatibility  of  Tables  lA  and  IB  is  only  general. 


Table  IB.  Near-field  anomaly,  decibels  (sea  state  2); 
ky  =  l0logKi;ky  =  \0  \ogKy 


Frequency 

(Hz) 

Sand 
kt  ky 

Mud 

ky  ky 

112 

6.3 

11.8 

6.3 

11.3 

446 

6.1 

10.3 

5.8 

9.1 

1120 

5.0 

6.7 

4.5 

5.8 

2820 

3.7 

4.3 

3.3 

3.8 
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D.  EfTective  Attenuation  Rate  —  Surface  Scattering  and  Bottom  Coupling 

The  next  consideration  to  be  discussed  is  a,,  the  effective  shallow- water  attenuation 
coefficient  in  decihels/skip  distance.  An  energy  conservation  rule  was  used  to  establish 
this  coefficient,  which  includes  the  additional  loss  due  to  the  coupling  of  energy  from  the 
wind-roughened  sea  surface  to  the  bottom. 

M-S  postulated  that  if  r,  is  the  surface  reflection  coefficient  then  (1-r,)  is  the 
surface-scattering  coefficient,  and  a,  =  -10  logr,  is  the  surface  loss  in  decibels/bounce. 

Sea  surface  scattering  has  three  effects  on  attenuation  rate  in  the  channel: 

(1)  It  smooths  the  depth  dependence  in  the  channel. 

(2)  It  modifies  the  near-field  anomaly. 

(3)  It  increases  the  rate  of  mode  stripping  so  that  the  range  to  single-mode 
transmission  (cylindrical  spreading)  is  shortened  considerably  to  8//, 
where  H  is  skip  distance. 

Figure  2  is  a  plot  of  sea  surface  loss  versus  frequency  x  wave-height  product  found 
from  analysis  of  the  AMOS  data  for  the  various  sea  states  and  corresponding  wind  condi¬ 
tions.  The  equations  in  Fig.  2  are  based  on  Ref.  9.  Figure  3  shows  the  same  plot  in 
metric  units*®  with  an  analytic  expression  for  wave  height** 

The  surface  components  r,  and  (1-r,)  must  undergo  different  interactions  with  the 
bottom  if  the  sea  surface  is  not  flat.  The  simplest  expression  satisfying  this  requirement 
and  the  Colossus  data  is 

where  r,  is  the  fraction  of  energy  transmitted  down  the  channel  when  a  bottom  event  is 
coupled  with  each  surface  reflection,  and  r^,  is  the  bottom  reflection  coefficient  The 
shallow-water  attenuation  coefficient,  in  decibels/bounce,  would  then  be 

a,  =  -10  logr,  . 

The  fact  that  must  be  multiplied  by  indicates  that  near-grazing  rays  suffer  a  bottom 
loss.  The  fact  that  the  scattered  rays  must  be  multiplied  by  means  that  the  angular 
dependence  of  the  bottom  reflection  loss  is  such  that,  on  the  average,  these  steeper  rays 
suffer  twice  the  loss  of  the  near-grazing  rays. 

The  shallow-water  attenuation  coefficient  values  for  sand  and  mud  bottoms  are 
given  as  a  function  of  sea  state  and  frequency  in  Table  n. 
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Figure  3.  The  dependence  of  acoustic  loss  at  the  sea  surface  (Ref.  10)  on  wind  speed. 


Table  II.  Shallow-water  attenuation  a,  (decibelsibounce). 


Sea 

Sute 

/ 

(kHz) 

0 

Sand 

Mud 

1 

Sand 

Mud 

2 

Sand 

Mud 

3 

Sand 

Mud 

4 

Sand 

Mud 

5 

Sand 

Mud 

0.1 

1.0 

1.3 

1.0 

1.3 

1.0 

1.3 

1.0 

1.3 

1.0 

1.3 

1.0 

1.3 

0.2 

1.3 

1.7 

1.3 

1.7 

1.3 

1.7 

1.3 

1.7 

1.3 

1.7 

1.4 

1.7 

0.4 

1.6 

2.2 

1.6 

2.2 

1.6 

2.2 

1.6 

2.2 

1.7 

2.4 

2.2 

3.0 

0.8 

1.8 

2.5 

1.8 

2.5 

1.9 

2.6 

2.2 

3.0 

2.4 

3.8 

2.9 

4.0 

1.0 

1.8 

2.7 

1.9 

2.7 

2.1 

2.9 

2.6 

3.7 

2.9 

4.1 

3.1 

4.3 

2.0 

2.0 

3.0 

2.4 

3.5 

3.1 

4.4 

3.3 

4.7 

3.5 

5.0 

3.7 

5.2 

4.0 

2.3 

3.6 

3.5 

5.2 

3.7 

5.5 

3.9 

5.8 

4.1 

6.2 

4.3 

6.4 

8.0 

3.6 

5.3 

4.3 

6.3 

4.5 

6.7 

4.7 

6.9 

5.0 

7.3 

5.1 

7.5 

10.0 

4.0 

5.9 

4.5 

6.8 

4.8 

7.2 

5.0 

7.5 

5.2 

7.8 

5.3 

8.0 
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m.  ROGERS*  ONBOARD  PREDICTION  MODEL  FOR  PROPAGATION 
LOSS 

A.  General 

Dr.  PJH.  Rogers  of  the  Naval  Research  Laboratory  prepared  a  report^  that  reviewed 
the  infonnation  base  for  shallow-water  acoustic -propagation  loss,  especially  in  the  In- 
quency  range  100  to  800  Hz.  One  of  his  chief  contentions  was  that  so  many  parameters 
("no  fewer  than  24")  were  required  to  determine  the  shallow-water  propagation  loss  that 
it  was  easy  to  explain  measured  results  but  hard  to  predict  them.  We  quote  his  complete 
abstract 

This  report  examines  the  state  oT  the  art  in  the  prediction  of  propagation  loss  in 
shallow  water  as  it  pertains  to  onboard  performance  prediction.  The  following  conclu¬ 
sions  are  drawn:  1.  For  simple  cases,  Le.,  homogeneous  liquid  bottom,  linear  sound- 
speed  gradient  no  surface  or  tottom  roughness,  a  simple  algebraic  model,  for  depth  aver¬ 
aged  propagation  loss  works  as  well  as  the  mote  complex  mode  nxxlel.  (The  model  is 
presented  in  the  report);  2.  The  uncertainty  in  bottom  parameters,  particularly  sound 
velocity  and  attenuation  makes  it  impossible  to  set  meaningful  bounds  on  fx-opagation 
loss  particularly  for  negative  gradients  or  slow  bottoms.  (Useful  predictions,  however, 
can  probably  be  nude  when  a  positive  gradient  is  present);  3.  Details  of  the  sound-speed 
profile  can  cause  significant  changes  in  propagation  loss,  therefore  even  if  bottoms  were 
well  characterized,  sophisticated  computer  models  would  be  required  to  predict  propaga¬ 
tion  loss;  4.  Virtually  all  propagation  loss  curves  can  be  described  to  within  a  fraction  of 
a  dB  [in  the  range  5  to  1(X)  km]  by  the  function  FL  *  B  +  15  logB  +AJt  +CJi^ 
with  the  C  coefficient  usually  zero.  Thus,  the  output  field  can  be  described  by  two  or,  at 
most  three  free  parameten.  Since  there  are  no  fewer  than  24  input  parameters  it  is  thus 
easy  to  explain  observed  propagation  loss  and  very  difficult  to  predict  it  Moreover,  it  is 
doubtful  that  propagation  loss  experiments  can  uniquely  define  bottom  parameters;  S. 

Certain  aspects  of  the  theory  remain  unverified  and/or  inadequately  treated.  These 
include:  (1)  surface  and  bottom  roughness.  (2)  shear  in  the  sediment,  (3)  substrate  rough¬ 
ness,  (4)  modal  coupling,  and  (S)  biological  scatterers;  6.  Grain  size  distribution  is  not 
an  adequate  predictor  of  acoustical  properties;  hence  currently  existing  sediment  charts 
are  of  little  or  no  value  in  performance  prediction;  and  7.  Many  input  parameters  arc 
very  poorly  known.  These  include;  a.  bottom  roughness,  b.  wave  height  spectrum,  c. 
sediment  shear-wave  speed,  d.  sediment  shear  attenuation,  e.  shear  and  sound-speed  and 
attenuation  gradients  in  the  sediment,  and  f.  distribution  and  effective  attenuation  of 
biologies.  In  most  cases,  the  theory  is  not  certain  enough  to  determine  the  uncertainty  in 
propagation  loss  caused  by  uncertainty  in  these  parameters. 

We  will  discuss  his  conclusions  as  applicable  to  the  M-S  model,  but  not  in  his  order. 

First  we  would  like  to  add  two  more  inputs  to  his  list  (Table  III)  of  "no  fewer  than 
24  inputs"  to  his  Universal  (Range  Independent)  Shallow  Water  Propagation  Loss  Model. 
The  overabundance  of  possibly  required  inputs  is  certainly  evident  from  Table  III.  Our 
question  is,  Can  we  selectively  reduce  this  number  to  a  manageable  few  which  are 
important  for  specific  categories  of  problems? 
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Table  HI.  Inputs  to  Universal  (Range  Independent)  Shallow  Water  Propagation  Loss  Model 


1.  Water  depth 

2.  Sound  speed  profile 

A.  Temperature 

B.  Salinity 

Water  column 

3.  Acoustic  attenuation  in  water 

4.  Internal  waves  and  tides 

5.  Sloping  bottom 

6.  Density  of  sediment 

7.  Sound  speed  in  sediment 

8.  Shear  speed  in  sediment 

9.  Acoustic  attenuation  in  sediment 

10.  Shear  attenuation  in  sediment 

11.  Sound  speed  gradient  in  sediment 

12.  Shear  speed  gradient  in  sediment 

13.  Attenuation  gradient  in  sediment 

14.  Density  gradient  in  sediment 

15.  Thickness  of  sediment  layer 

X  number  of  layers 

16.  Sound  speed  in  basement 

17.  Shear  speed  in  basement 

18.  Density  of  basement 

19.  Acoustic  attenuation  in  basement 

Basement 

20.  Shear  attenuation  in  basement 

21.  Surface  roughness 

22.  Bottom  roughness 

23.  Subbottom  roughness 

24.  Entrained  gas  bubbles 

25.  Fish  and  other  biological  scatterere 

26.  Wind  vector 

Scattering 

His  report  recognizes  that  "Useful  predictions,  however,  can  probably  be  made 
when  a  positive  gradient  is  present."  This  is  precisely  the  situation  depicted  by  the  M-S 
model. 

We  note  that  his  model  is  a  semiempirical  one  based  on  computed  normal-mode 
solutions.  In  the  report,  he  covers  only  profiles  in  which  the  sound  speed  decreases  mono- 
tonically  with  depth  and  the  average  gradient  is  0.2  s~*.  In  addition,  he  covers  only  situa¬ 
tions  without  sea-surface  scattering. 
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Rogers  states  that  virtually  all  computed  (depth  averaged)  propagation  loss  curves 
can  be  fit  to  within  a  fraction  of  a  decibel  to  the  following  simple  function  (for  most 
cases  C=0): 

PL -ISlogR +AR +B +CR^  (12) 

where  PL  is  the  propagation  loss  in  decibels  and  R  is  the  range  in  kilometers  (5  to 
100  km). 

Since  only  two  free  empirical  algebraic  parameters,  A  and  B,  are  needed  to  fit  his 
equation  between  5  and  100  km,  we  will  relate  these  free  parameters  to  his  theoretical 
expression  and  the  M-S  theoretical  base. 

Rogers’  theoretical  model  for  propagation  loss  in  water  with  a  strictly  negative 
sound-speed  profile  leads  to  his  Eq.  (4)  (Eq.  (13)  here].  (Weston^  has  develored  an 
analogous  model  for  the  isovelocity  case.) 

PL=\SlosR  +5log(H^)  +  ^-^-lAZ  +  a^R  dB,  (13) 

AH 


where  15  logf?  {R  is  in  meters)  is  the  spreading-loss  term  for  the  mode- stripping  region, 
H  is  the  water  depth  in  meters,  ^  is  bottom  loss  in  decibels/radian  and  comes  from  the 
theoretical  expression  for  the  Rayleigh  reflection  coefficient  for  a  two-fluid  lossy  inter¬ 
face  at  small  grazing  angles,  and  Oy,  is  the  absorption  coefficient  of  seawater  in 
dccibels/kilomcter  and  is  common  to  all  models.  is  the  larger  of  0^  or  0^ ,  where  Bg 
is  the  maximum  grazing  angle  for  a  skip  distance.  The  definitions  of  Bg  and  6^  follow: 

0,  =yl(ZHgf,VC^  .  (14) 


where  C„  is  the  value  of  the  sound  speed  (maximum)  at  the  surface  of  the  water  and  gfj 
is  the  magnitude  of  the  negative  sound-speed  gradient,  in  meters/second/meter,  or 
seconds"^;  0^  is  in  radians.  0^.  is  the  effective  plane- wave  angle  for  the  lowest  propagat¬ 
ing  mode. 


(15) 


where  /  is  the  frequency  in  hertz,  and  0^  is  in  radians. 

For  most  cases  of  interest,  Eller’s  expression  was  used  by  Rogers^  to  obtain  the 
reflection  coefficient 


OAnMpNQk^ 


dB/rad , 


(16) 
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where 

Nq  =  CJCj  and  A/o  =  P,/pw 

=  maximum  (sea  surface)  sound  speed  in  the  water 
Cj  =  sound  speed  in  the  sediment 
s=  density  of  the  water 
p,  =  density  of  the  sediment 

kp  =  sediment  attenuation  coefficient,  decibels/meter/kilohertz. 


B.  Frequency  Dependence 

Rogers  claims  that  for  negative  sound-speed  gradients  the  effective  attenuation 
coefficient  A  is  independent  of  frequency  over  a  wide  range  of  frequencies.  The  assump¬ 
tion  in  this  statement  is  that  kp  is  constant  over  this  frequency  range,  corresponding  to 
the  Hamilton  loss  factor.  This  assumption  has  been  accepted  widely,  but  not  uncondition¬ 
ally  or  universally.  This  is  an  important  issue,  because  the  M-S  model,  which  is  based  on 
measurement  data,  shows  a  distinct  frequency  dependence  due  to  both  sea  surface 
scattering  loss  and  bottom  loss.  In  fact,  kp  seems  to  vary  with  frequency  as  a  higher 
power  than  1  and  may  even  be  as  much  as  2.  Ingenito^^  found  kp  to  vary  as  /  in  the 
Gulf  of  Mexico  off  Panama  City.  In  other  work,*^  the  measured  frequency  dependence 
of  the  attenuation  term  has  been  used  as  a  clue  to  the  nature  of  the  attenuation  process 
controlling  the  propagation  at  a  particular  location  and  time. 

C.  Comparison  of  Rogers’  and  M-S  Equations 

Consider  Rogers’  empirical  Eq.  (17)  for  propagation  loss  in  the  mode-stripping 
region  with  the  very  small  term  CR  ^  dropped: 

PL  =  15  log/? +/y? +B  dB.  (17) 

Rogers  estimates  A  and  B  by  solving  the  wave  equation  by  normal  modes  for  a  specific 
set  of  conditions.  He  then  fits  the  computed  points  by  the  least-squares  method. 

If  we  write  Rogers’  solution  [Eq.  (13)]  for  propagation  loss  under  strictly  negative- 
gradient  conditions  and  also  the  M-S  Eq.  (2)  for  the  mode- stripping  region,  we  find  that 
both  equations  can  be  related  to  A  and  B  of  Eq.  (17). 
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Rogers  [Eq.  (13)]: 


Pi,  =  151ogR  +51og(HP)  +  -^« -7.18  +  a„R 

4// 


M-S  [Eq.  (2)1: 

A/  =  15  log/?  +  5  log//  +  a, 


r  ^ 


-kl  +  60  +  a/? 


dB 


dB 


for//  <.R  S8//. 


In  these  equations,  a,,,  (dB/m)  and  a  (dB/kyd)  are  the  same  absorption  coefficient.  Note 
also  that  in  Rogers’  equation  H  stands  for  water  depth,  while  M-S  use  H  for  skip  dis¬ 
tance. 


Equation  (13)  holds  from  the  end  of  the  direct  radiation  zone  to  range  R  jq,  which  is 
where  the  10  log/?  zone  begins. 


R 


3/2 


pc2 


(18) 


Here,  H  is  the  water  depth  and  P  is  given  by  Eq.  (16).  Note  that  R  jo  is  now  a  function  of 
frequency. 

The  corresponding  coefficients  in  Eqs.  (13)  and  (2)  have  similar  meanings.  There  is 
an  attenuation-rate  term  multiplying  the  range  ai'.d  a  constant  term  representing  the  refer¬ 
ence  level  at  which  energy  begins  to  feed  into  the  channel  from  the  direct  radiation  zone. 
A  major  difference  concerns  the  extent  of  the  mode-stripping  region.  Rogers  extends  it 
from  5  km  to  100  km  instead  of  using  Eq.  (18)  where  bottom  losses  represented  by 
Eller’s  expression  for  P  may  cause  it  to  shorten.  The  use  of  P  is  an  approximation  that 
holds  for  small  grazing  angles  only. 

The  M-S  model  uses  a  comparatively  short  mode-stripping  region,  H  to  8// ,  based 
on  measurements  off  the  Atlantic  Coast.  The  limiting  value  of  %H  is  mainly  due  to  the 
low-frequency  data,  which  showed  relatively  large  losses  at  frequencies  of  about  lOO  Hz 
to  500  Hz.  These  large  losses,  as  well  as  differences  in  bottom-type  classification,  are 
fundamental  to  the  differences  in  the  values  predicted  by  the  M-S  and  Rogers  models. 
Because  sea-surface  scattering  is  small  at  low  frequencies,  the  larger  losses  must  be  due 
to  interaction  with  the  bottom. 

The  properties  of  the  subbottom  also  become  important  at  low  frequencies.’^  The 
M-S  model  assumes  that  the  bottom  is  homogeneous.  The  bottom  can,  however,  contain 
posuive  sound-speed  gradients  which  can  affect  propagation  drastically  by  returning 
energy  to  the  water.  For  example,  the  penetration  depth  into  clayey  silt  at  200  Hz  for  the 
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particular  case  under  study  can  be  as  much  as  23  m;  for  silt,  the  penetration  depth  is  4  m, 
and  for  fine  sand,  3  m.  To  depths  of  1  to  2  m  beneath  the  seabed  surface,  the  sound-speed 
gradient  is  about  +1.0  s”^  for  silts  and  clays  and  about  +15  s"*  for  fine  sand.*^ 
Substructure  layers  can  also  affect  energy  returned  to  the  water,  or  the  apparent  bottom 
loss.  Long  bottom  paths  can  feed  energy  trapped  in  the  bottom  back  into  the  shallow- 
water  channel  anywhere  along  the  path  depending  on  the  lateral  inhomogeneities  and 
slopes  of  sedimentary  layers.  In  South  China  Sea  areas,  for  example,  workers  have  found 
the  M-S  low  frequency  predictions  for  the  Atlantic  Coast  give  losses  that  are  too  high.^^ 

Rogers  is  able  to  test  the  sensitivity  of  his  model’s  predictions  to  changes  in  param¬ 
eters  such  as  bottom  roughness  or  the  positive  and  negative  sound-speed  gradi'mts  in  the 
bottom.  He  finds  that  there  can  be  very  large  changes  in  the  resulting  curves  of  propaga¬ 
tion  loss  vs  range.  Thus  there  is  an  important  requirement  to  measure  and  know  the  pro¬ 
perties  of  the  bottom  in  the  particular  geographic  area  of  interest. 

Rogers’  overall  point  is  that,  for  propagation  loss  in  shallow  water,  there  are  many 
parameters  that  can  be  used  to  develop  equations  for  A ,  the  attenuation  rate,  and  B ,  the 
constant  term  corresponding  to  contributions  by  the  near-field  anomaly.  To  test  these 
parameters,  he  chooses  to  use  a  theoretical  and  computational  method  based  on  normal¬ 
mode  solutions  to  the  wave  equation.  This  approach  may  be  appropriate  for  sensitivity 
analyses,  but  any  absolute  use  must  be  backed  by  measurements  of  the  parameters  that  he 
considers. 

Ttie  M-S  approach  was  quite  specific  in  developing  a  semiempirical  model  based  on 
observed  physical  data.  Although  originaUy  designed  to  handle  average  sound-speed 
structures  with  variability  in  bottom  type,  wind  speed,  and  layer  depth,  the  M-S  model 
can  now  be  adapted  to  specific  negative-gradient  profiles  like  those  covered  by  Rogers’ 
model. 

On  the  other  hand,  normal-mode  solutions  have  been  obtained  with  boundary 
scattering^^  and  can  be  adapted  to  Rogers’  method.  Unless  tire  solutions  are  tested  by  a 
measurement  program  and  error  bounds  estabUshed  such  as  in  the  M-S  model,  the  model 
cannot  be  considered  validated.  Next  we  compare  the  results  of  the  two  models  and 
show  that  they  are  significantly  different 

D.  Comparison  of  Predicted  Propagation  Losses 

One  of  the  purposes  of  the  present  study  is  to  compare  the  results  of  the  M-S 
shallow-water  prediction  method  developed  in  1962  with  the  onboard  prediction  method 
devised  20  years  later  (1981)  at  NRL^  using  empirical  fits  of  normal-mode  solutions  by 
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frequency  and  bottom  type.  In  Table  IV,  predictions  of  propagation  loss  vs  range  are 
listed  for  Rogers’  fine  sand  and  clayey  silt  at  200  Hz.  Alongside  are  placed  the  M-S 
predictions  for  the  Colossus  "sand"  and  "mud."  These  latter  classes  are  more  general,  and 
there  is  no  reason  to  believe  they  correspond  to  the  classes  for  wh’  Th  Rogers  calculated 
predictions. 


Table  TV.  Comparison  of  Rogers’  and  M-S  predictions  for  propagation 
loss  (decibels)  ai  200  Hz. 


Range 

(km) 

Fine  Sand 
Rogers 

Sand 

M-S 

A 

Clayey  Silt 
Rogers 

Mud 

M-S 

A 

5 

60 

69 

+9 

68 

69 

+1 

10 

66 

74 

+8 

76 

75 

•  1 

20 

72 

81 

+9 

89 

82 

-7 

3C 

77 

86 

+9 

99 

88 

-11 

40 

81 

90 

+9 

108 

93 

-15 

50 

84 

94 

+10 

117 

97 

-20 

60 

87 

97 

+10 

126 

101 

-25 

70 

90 

100 

+10 

135 

104 

-31 

80 

93 

103 

+10 

143 

108 

-35 

90 

96 

105 

+9 

152 

111 

-41 

100 

99 

108 

+9 

160 

115 

■45 

Table  IV  shows  that,  for  a  fine*sand  bottom,  the  M*S  "sand"  model  predicts  a 
greater  loss  than  Rogers’  fine-sand  model  by  about  9  or  10  dB  at  all  ranges.  On  the  other 
hand,  the  M-S  mud  model  predicts  a  smaller  loss  than  Roger’s  clayey-silt  model,  with  the 
difference  increasing  with  range.  The  differences  could  correspond  to  differences  in  clas¬ 
sifying  the  bottom  or  they  could  be  due  to  the  prediction  technique  or  both.  The  differ¬ 
ences  A  far  exceed  the  M-S  probable  errors  shown  in  Table  V.  Where  the  difference  is 


Table  V.  Probable  error  of  propagation  loss  (decibels) 
(semi-interquartile  range). 


Range 

(kyd) 

112 

Frequency  (Hz) 
446  1120 

2820 

3 

2 

4 

4 

4 

9 

2 

4 

5 

6 

30 

4 

9 

11 

11 

60 

5 

9 

11 

12 

90 

6 

9 

11 

12 
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constant  like  that  for  the  sands,  it  could  be  due  simply  to  the  value  used  for  the  bottom 
reflection  coefficient  or  to  an  error  in  the  near-field  anomaly  estimate.  Where  the  differ¬ 
ences  increase  with  range  as  for  the  muds,  they  are  probably  due  to  the  estimate  of  the 
attenuation  rate.  There  are  nine  classes  in  Hamilton’s  system,  but  Rogers  does  not  tabu¬ 
late  prediction  fits  for  other  classes  of  sand  and  silt 

Table  VI  compares  the  environmental  and  acoustic  parameters  used  by  the  two 
models.  The  most  important  differences  are  the  skip  distance  and  the  reflection  loss  per 
bounce.  Note  that  for  the  sands,  the  attenuation  rate  is  about  the  same  in 
decibels/kilometer  because  it  is  equal  to  the  quotient  of  loss  per  bounce  by  skip  distance. 
This  is  not  true  for  the  muds  or  silts.  Both  the  skip  distance  and  the  bottom  reflection 
loss  per  bounce  depend  on  the  magnitude  of  the  negative  gradient  from  the  surface  to  the 
bottom. 

The  question  to  be  answered  is.  Can  we  extend  the  M-S  model  to  predict  the  same 
propagation  loss  as  Rogers’  normal  mode  procedure  if  the  same  shallow-water  charac¬ 
teristics  are  used?  We  will  show  that  the  answer  is,  yes. 

Table  VI.  Environmental  acoustic  parameters  for  the  two  models 
(f  =200  Hz). 


Parameter 

Rogers 

Model 

M-S 

Model 

I.  Environment 

Layer  depth,  m 

0 

0 

Water  d^th,  m 

100 

100 

Negative  gradient,  s"' 

-0.20 

-0.035 

Sea  state 

C 

0 

Skip  distance,  km 

2.45 

5.86 

Skip  ray  grazing  angle,  deg 

9.36 

3.91 

II.  Reflection  Loss 

Fine  sand: 

dB /bounce 

0.48‘ 

1.3 

dB/km' 

0.20 

0.22 

/?io.  km 

100 

47 

Clayey  silt: 

dB /bounce 

3.46‘ 

1.7 

dB/km' 

1.41 

0.29 

/?,o,  km 

100 

47 

'Water  absorption  coefficient  must  be  added  to  this  term 
It  is  4.63  X  10"’  dB/km. 


‘Calculated  from  Morsc-Mackenzie. 
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IV.  EXTENSION  OF  THE  M-S  MODEL 

A.  Arbitrary  Negative  Sound>Speed  Gradients 

The  M-S  model  shown  in  Table  VI  uses  the  same  value,  -0.035  s"^  for  all  negative 
gradients.  We  have  developed  an  extension  to  the  M-S  model  for  handling  arbitrary 
negative  gradients  like  the  -0.2  s~^  that  Rogers  uses.  Changing  the  gradient  will  affect 
both  the  reflection  loss  per  bounce  and  the  near-field  anomaly  and  will  result  in  values 
comparable  to  Rogers’  since  the  skip  distance  and  grazing  angle  will  also  be  comparable. 

Near -Field  Anomaly 

The  general  expression  for  the  near-field  anomaly  ki  y  is  obtained  by  combining 
Eqs.  (6)  and  (9). 

.  ..s.  ^  L  .  Tfc-l-r, -l-riXl-r^r^)  ] 

^L,U  ~  ~  ^  - T, - T ^  ^  • 

U  -'’s  ’'b) 


r,,rb  S  1  . 


The  use  of  Ki^  y  instead  of  Ky  or  Ky  is  one  of  the  keys  to  the  extended  M-S  model. 
In  particular,  we  must  determine  the  number  n  of  the  surface-reflected  bottom-reflected 
(SRBR)  ray  system  that  applies  in  a  particular  situation.  For  a  strictly  negative  sound- 
speed  gradient,  we  proceed  as  follows  (see  Fig.  4). 


SRBR  ray  • 


Critical 
/0A  angle 


D 

(depth) 


Maximum 
^  skip  ray 


ical 
angle  range) 


Order  n  =  Integer  (H/R  ) 


(skip  distance) 


Figure  4.  Determination  of  number  n  of  SRBR  rays  in  direct  radiation  zone. 
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First  we  state  that  the  critical  angle,  i.e.,  the  cutoff  angle  for  accepting  energy  into 
the  channel,  is  given  by 

Gj.  =  arc  cos  Nq  ; 

Nq  =  ci/c2,  where  c  i  is  the  sound  speed  in  the  water  at  the  bottom  and  c  2  is  the  sound 
speed  in  the  bottom  sediment 

In  this  situation  of  negative  gradients  and  low  sea  states,  surface  reflection  =  1. 
Thus  more  rays  at  higher  angles  up  to  the  critical  angle  contribute  to  t’  .e  field.  If  is  the 
range  between  bottom  contacts  of  the  ray  (not  refracted)  traveling  at  the  critical  angle 
with  the  bottom,  then,  from  simple  geometry, 

Rc -  2D  cot  0^  . 

The  number  of  bottom-reflected  rays  that  contribute  to  the  near-field  anomaly  is 
integer  n  ^HIR^,  where  H  is  the  skip  distance  (see  Section  II.B).  In  Table  Vn,  we 
show  how  to  obtain  n  for  water  with  a  depth  of  100  m  and  a  negative  sound-speed  gra¬ 
dient  of  -0.2  s”\  for  two  types  of  bottom  sediments  —  fine  sand  and  clayey  silt. 

Again,  we  repeat  that  the  determination  of  the  near-field  anomaly  is  not  trivial  since 
the  same  bottom  reflection  coefficient  r^,  must  be  applicable  to  all  three  propagation 
2X>ncs  —  near,  intermediate,  and  far. 

To  calculate  the  extended  near- field  anomaly,  ki^  y,  for  the  two  sediments,  we  need 
to  know  their  reflection  coefficients,  or  bottom  losses  per  bounce. 

Analysis  of  the  Colossus  data  yielded  the  bottom  attenuation  factors  vs  frequency 
for  sand  and  mud  shown  in  Fig.  5 ;  these  values  are  for  a  sea  state  0  and  the  shallow 
waters  off  the  East  Coast  of  the  United  States.  Using  the  Morse-Mackenaic**  equations 
for  lossy  sediments  modeled  as  fluids,  we  will  now  study  the  relation  between  these 
empirical  values  and  theoretical  values. 

Toole  VII.  Number  n  qfSRBR  rays  in  direct  radiation  zone.^ 


Parameter/Sediment 

Fine  Sand 

Clayey  Silt 

Skip  distance,  m 

2530 

2530 

Critical  angle,  O^,  deg. 

29.12 

9.45 

Critical-angle  ray  range,  Rc,  m 

359 

1191 

Order  of  ray  bounces,  n 

7 

2 

“Water  depth  is  100  m  and  sound  speed  gradient  is  -0.2  s"'. 
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where 

6  =:  grazing  angle,  degrees 
/  =  frequency,  hertz 

ag  =  attenuation  in  the  bottom,  decibels/meter 
vg  =  phase  velocity  in  the  bottom,  =  C2 


As  an  example,  using  the  bottom  properties  for  fine  sand  and  clayey  silt  as  defined 
by  Hamilton*’  from  Table  VIII,  we  obtain  the  curves  shown  in  Fig.  6  for  bottom  loss 
versus  grazing  angle  at  200  Hz. 

Table  VIII.  Bottom  properties  (Hamilton}’^}. 


Parameter 

Fine  Sand 

Clayey  Silt 

Seawater 

Sound  speed,  m/s 

1749 

1549 

1528 

Density,  g/m^ 

1.941 

1.488 

1.027 

Hamilton  attenuation 
factor,  dB/m/kHz 

0.52 

0.095 

... 

Critical  grazing  angle, 

A»rt 

29.12 

9.45 

We  have  pointed  out  in  discussing  the  near-field  anomaly  that  r^,  turns  out  to  be  the 
reflection  coefficient  of  rays  up  to  the  critical  angle.  This  is  appropriate  for  a  scattering 
or  reflection  process  in  a  channel  since  this  is  the  cutoff  angle  for  accepting  energy  into 
the  channel. 

E.L.  Hamilton*’  has  postulated  from  many  measurements  that  the  loss  for  compres- 
sional  waves  propagating  in  a  specific  type  of  sediment  has  a  fixed  attenuation  constant, 
kp ,  in  decibels/meter/kilohertz.  If  this  were  true,  the  Morse-Mackenzie  expressions 
would  show  that  the  bottom  loss  curve  is  independent  of  frequency.  Figure  5  shows  that 
this  last  condition  is  not  true  for  the  Colossus  data  for  sand  or  for  mud. 
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Figure  6.  Bottom  loss  (decibels)  vs  grazing  angle  (degrees)  (Morse-Mackenzie). 


On  the  other  hand,  if  we  take  kp  for  fine  sand  and  also  clayey  silt  to  vary  as  the  first 
power  of  frequency  up  to  1  kHz  and  calculate  the  loss  per  bounce,  then  the  corresponding 
angle  of  incidence  is  approximately  the  critical  angle  for  each  sediment.  Above  1  kHz, 
taking  kp  equal  to  a  constant  also  yields  a  bottom  loss  corresponding  to  the  critical  angle 
(see  Table  IX).  This  is  reasonable,  since  the  critical  angle  is  the  cutoff  angle  for 
transmission  down  the  channel,  and  the  scattering  is  quite  diffuse.  Often,  Hamilton’s  pos¬ 
tulate  is  assumed  for  convenience,  as  Rogers  does  in  his  model.  However,  some  other 
workers  have  also  found  a  behavior  similar  to  that  of  the  Colossus  data  for  the  frequency 
dependence  of  kp 

TcUfle  DC.  Empirical  bottom  loss  per  bounce  (MS)  and  corresponding  grazing  angle. 


Frequency 

(kHz) 

Sand 

(dB/bounce) 

e; 

(deg) 

Mud 

(dB/bounce) 

e; 

(deg) 

0.1 

1.0 

29.4 

1.3 

9.47 

0.2 

1.3 

29.6 

1.7 

9.46 

0.4 

1.6 

29.5 

2.2 

9.40 

0.8 

1.8 

29.3 

2.5 

9.01 

1.0 

1.8 

28.6 

2.7 

8.82 

2.0 

2.0 

29.0 

3.0 

9.06 

4.0 

2.3 

29.5 

3.6 

9.44 

Median  -  29.3 

Median  -  9.40 

•0,  is  calculated  from  Morse-Mackenzie  for  •«  /  when  /  s  1  kHz; 
kp  is  constant  when  /  a  1  kHz  (see  p.  M  for  definition  of  kp). 
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Skip  Distance 


A  new  expression  for  the  skip  distance  can  be  derived  based  on  the  actual  gradient 
in  the  negative-gradient  layer  rather  than  on  the  average  for  the  Colossus  measurements. 
At  this  point,  we  are  concentrating  on  the  strictly  negative-gradient  profile.  (The  formula 
for  a  bilinear  gradient  is  provided  in  Section  IV.B.)  For  strictly  negative  sound-speed  gra¬ 
dients,  there  are  two  useful  expressions;  (1)  maximum  skip  angle,  6^  and  (2)  skip 
distance,  H. 


e 


g  max 


^8n 


*=1 


1/2 


radians 


and 


H  =2 


2Dc 


m 


gN 


1/2 


where  D  is  the  water  depth,  gfj  is  the  magnitude  of  the  negative  gradient,  c  j  is  the  sound 
speed  of  the  water  at  the  bottom,  and  is  the  maximum  sound  speed  (at  the  surface). 


When  only  negative  gradients  exist,  the  sea  state  must  be  close  to  0,  for  otherwise 
near-surface  mixing  by  the  wind  would  begin  to  form  isothermal  layers.  Under  this  con¬ 
dition,  the  sea  surface  reflection  coefficient  r,  is  1,  and  r,  =  r^. 

As  an  illustration,  we  show  in  Table  X  the  near-field  anomaly  ki^  y  for  two  of 
Hamilton's  sediment  classes  for  shallow  water  of  gf^  =  0.2  s"^  D  *  1(X)  m,  //  =  2.45  km, 
and  0,  =  9.356'’. 


Table  X.  Determination  of  near-field  anomaly  for  two  Hamilton  sedment  classes 
at  200  Hz. 


Property/Sediment 

Fine  Sand 

Clayey  Silt 

Reflection  coefficient,  r^, 
at  skip  angle 

0.90 

0.45 

Bottom  loss,  dB/bounce 

0.48 

3.46 

Critical  Siigle,  6,,,  deg. 

29.12 

9.45 

Critical  cycle  range,  Rc,  m 

359 

1191 

Order  of  ray  bounces,  n 

7 

2 

Near-field  anomaly,  dB 

13.1 

7.9 

gs  =  0.2  s"';  D  =  100  m;  H  (skip) 

=  2.45  km;  O, 

=  9.356“ 
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At  this  point,  it  is  useful  to  discuss  the  shallow-water  propagation-loss  measure¬ 
ments  of  Cole  and  Podcszwa^*  south  of  Long  Island  in  30  fathoms  (54.9  m)  of  water  over 
a  sand  bottom  at  3.5  kHz  under  strictly  negative-gradient  conditions.  The  significance  of 
this  experiment  is  that  it  was  devoted  to  exploration  of  the  details  of  the  standard  M-S 
model  under  such  conditions.  In  this  case,  they  found  that  M-S  got  the  right  answer  for 
the  wrong  reason  —  that  is,  use  of  the  average  Colossus  negative  gradient  gave  a 
predicted  result  within  the  probable  error  bounds  (see  Table  V),  but  the  Colossus  gradient 
did  not  correspond  to  the  actual  gradient  The  M-S  model  predicted  values  of  79  dB  vs 
their  measured  values  of  80  dB  at  10  kyd,  and  92  dB  vs  98  dB  at  20  kyd;  however,  Cole 
and  Podeszwa  found  a  skip  distance  of  1.20  kyd  and  a  mean  boundary  loss  of 

I. 3  dB/bounce  compared  with  the  skip  distance  of  4.75  kyd  and  the  mean  boundary  loss 
of  3.55  dB/bounce  predicted  by  the  M-S  model.  They  also  found  a  grazing  angle  of 

I I. 908°  at  their  skip  distance. 

To  show  the  versatility  of  the  basic  formulation  of  the  M-S  propagation-loss  equa¬ 
tions,  we  have  used  the  specific  environmental  parameters  found  in  the  Cole-Podcszwa 
experiment  and  have  calculated  predicted  propagation  loss  values  to  compare  with  their 
measured  values.  We  find  as  they  do  a  skip  distance  of  1.20  kyd  and  a  skip  angle  of 
1 1.908°  based  on  ray  traces  using  their  typical  sound-speed  profile.  From  their  values  of 
sand  density,  sound  speed,  and  a  sediment  attenuation  factor  of  0.422  dB/m/kHz  at 
3.5  kHz,  we  then  calculate  a  reflection  coefficient  of  0.75  (or  -1.25  dB/bounce).  Using  the 
reflection  coefficient  of  0.75,  we  calculate  ki  u  -  6.28  dB.  Then,  using  M-S  Eqs.  (2)  and 
(3),  we  get  a  propagation  loss  at  10  kyd  of  80  dB  vs  their  80  dB,  and  at  20  kyd  we  get 
96  d3  vs  their  98  dB. 

Comparing  the  Extended  M-S  Model  for  Arbitrary  Negative  Gradients  with  Rogers’ 
Model 

The  three  propagation-loss  zones  for  the  M-S  model  extended  to  include  arbitrary 
negative  gradients  are  as  follows; 

A/,dB  =  201og/? +a/? +60-k/,  (/  for/?  S//  (25) 

p 

A/,  dB  =  15  log/? +a/? +a,  (-— -1)  +  5  log// +  60 {/  for//S/?^jo  (26) 

A/,  dB  =  10  log/?  a/?  a,  (*^“1)  +  5  log//  +60  +  5  log/?  jq  -  (27) 

H 

for/?  S/?  10 , 
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where 

/?  10  =  range  to  the  start  of  the  far  region.  R  lo  is  obtained  from  Eq.  (18). 

a,  -  Of,  for  strictly  negative  sound-speed  gradients  from  surface  to  bottom. 

We  will  now  use  these  new  expressions  to  examine  the  relationships  between  the 
extended  M-S  and  Rogers’  fits  for  two  different  sediments,  fine  sand  and  clayey  silt. 

A  problem  that  we  encounter  immediately  is  that  Rogers  apparently  did  not  use  the 
Morse-Mackenzie  expression  for  the  reflection  coefficient  but  rather  Eller’s  expression, 
Eq.  (16).  Eller’s  expression  holds  only  for  very  small  grazing  angles.  For  clayey  silt,  it 
yields  a  reflection  loss  of  1.99dB/bouncc  instead  of  3.46  dB/bounce.  For  fine  sand,  it 
yields  a  reflection  loss  of  0.58  dB/bounce  instead  of  0.48  dB/bounce.  Using  these  Eller 
reflection  losses  in  Eqs.  (25),  (26),  and  (27),  we  obtain  the  results  shown  in  Table  XI. 
Thus  the  extended  Marsh-Schulkin  model  will  give  practicaUy  the  same  results  as 
Rogers’  if  the  same  inputs  are  used. 


Table  XI.  Comparison  of  Rogers’  and  extended  M-S  predictions  for  propagation 
loss  (decibels)  at  200  Hz. 


Range 

(km) 

Fine  Sand  ^ 

Rogers  Extended  M-S 

Clayey  Silt  ^ 

Rogers  Extended  M-S 

5 

60 

60 

0 

68 

66 

-2 

10 

66 

66 

0 

76 

75 

-1 

20 

72 

73 

*1 

89 

87 

-2 

30 

77 

78 

+1 

99 

98 

-1 

40 

81 

82 

+1 

108 

108 

0 

50 

84 

86 

+2 

117 

117 

0 

60 

87 

89 

+2 

126 

126 

0 

70 

90 

93 

+3 

135 

135 

0 

80 

93 

96 

+3 

143 

143 

0 

90 

96 

99 

+3 

152 

152 

0 

100 

99 

102 

+3 

160 

160 

0 

B.  Bilinear  Sound>Speed  Gradients 

The  extension  of  the  M-S  propagation-loss  model  is  based  on  the  three-zone 
spreading-loss  system  for  a  bilinear  positive-over-negative  sound-speed  gradient.  The 
original  M-S  model  was  a  three-zone  spreading-loss  system  derived  from  measured 
acoustic  data.  By  using  the  Hamilton  sediment  classes  for  estimating  the  acoustic  pro¬ 
perties  of  the  bottom  from  1(X)  Hz  to  5(X)0Hz,  we  obtain  the  reflection  loss  per  bou-.ice 
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from  one  of  the  loss-vs-grazing*angIe  expressions  developed  by  Morse-Mackenzie,^® 
Brekhovskikh,^  Weston,^  or  Eller.  The  Schulkin-Marsh  surface>scattering  loss  is  also 
used  for  a  given  wind  speed  or  sea  state.  The  water  absorption  coefficient  can  be  that  of 
Thorp^  or  Francois  and  Garrison.^ 


i 


The  first  rone  is  the  direct  path  zone  (20  log  R),  the  extent  of  which  is  defined  as  the 
skip  distance.  The  skip  distance  can  be  obtained  by  ray  tracing  or  from  the  formula 
given  below  for  bilinear  gradients: 

where 

=  maximum  sound  speed 
gp  =  positive  gradient 
gf^  =  magnitude  of  the  negative  gradient 
H  -  skip  distance 
D  =  water  depth 
L  -  layer  depth. 


jL+m=LL 

ip  8n 
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(28) 


The  near-field  anomaly  can  be  taken  as  ky  or  ,  given  by  Eq.  (9). 

The  second  zone  is  the  mode-stripping  zone  (IS  log  R)  in  which  the  shallow- water 
attenuation  coefficient,  a, ,  in  decibcls/bounce,  depends  on  the  grazing  angle  of  the  bot¬ 
tom  skip  ray  and  the  sea  state/wind  speed.  The  attenuation  coefficient  multiplied  by  the 
number  of  bounces  up  to  a  given  range  yields  the  attenuation  due  to  the  bottom  in  the 
second  and  third  zones. 
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The  third  rone  is  the  zone  of  cylindrical  spreading  (10  log  R)  in  which  there  is  only 
mode  of  propagation.  Instead  of  using  the  frequency-independent  empirical  value  of 
as  the  transition  range  between  zone  2  and  zone  3,  we  use  the  frequency-dependent 
given  in  Eq.  (18). 
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V.  ASSESSMENT  OF  ROGERS*  AND  M-S  MODELS 

Rogers’  method  of  fitting  two  (or  more)  coefficients  to  depth-averaged  normal¬ 
mode  solutions  of  the  wave  equation  for  acoustic  propagation  in  shallow  water  has  the 
advantage  of  providing  a  sensitivity  analysis  for  comparing  the  effects  of  adding  or 
changing  parameters  of  the  model.  In  addition,  the  equations  furnish  convenient  parame¬ 
ters  which  could  provide  reference  solutions  for  a  measurement  program. 

While  the  Rogers  models  are  computer-driven,  the  M-S  model  is  data-driven.  The 
latter  model  was  based  on  analysis  of  data  gathered  during  extensive  measurements  of 
propagation  loss  and  environmental  parameters  off  the  Atlantic  Coast,  and  basically 
applies  to  average  conditions  in  this  region.  To  use  it  anywhere  else  may  furnish  prelimi¬ 
nary  loss  estimates  but  is  not  strictly  justified.  The  main  reason  for  this  is  that  the  loss 
values  at  long  range  were  determined  mainly  by  the  low-frequency  properties  of  the  bot¬ 
tom  and  subbottom.  These  are  not  well  known  and  differ  from  one  location  to  another. 
Another  limiting  factor  to  its  application  is  that  it  concentrates  on  an  average  positive- 
gradient  layer  overlying  an  average  negative-gradient  layer  when  sea  surface  scattering 
might  be  present  When  sea  state  0  or  1  conditions  prevail  and  a  negative  gradient  is 
present,  it  only  gives  a  single  answer  for  all  conditions.  However,  we  have  shown  that  the 
M-S  method  can  be  extended  to  such  cases. 

Finally,  being  tied  to  a  measurement  data  base  is  an  important  advantage  for  the 
M-S  model  since  error  bounds  were  established  for  applicable  conditions. 
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VI.  CONCLUSIONS 

1.  The  M-S  shallow-water  transmission-loss  model  has  been  reviewed  and 
extended  in  the  light  of  20  years  additional  work  in  the  field.  It  appears  to  be  as  useful 
today  as  it  was  when  first  formulated.  Based  on  data  gathered  during  the  Colossus  pro¬ 
gram,  the  model  treats  surface-scattering  loss  coupled  to  the  bottom  for  a  water  column 
with  a  positive-gradient  layer  overlying  a  negative-gradient  layer  for  two  classes  of  bot¬ 
toms,  sand  and  mud. 

2.  In  terms  of  Hamilton’s  bottom  classes  and  associated  acoustic  properties,  the 
Colossus  sand  is  not  equivalent  to  his  "fine  sand"  and  the  mud  is  not  equivalent  to  his 
"clayey  silt."  Information  on  subbottom  structure  and  geoacoustic  properties  is  much 
needed.  Hamilton’s  postulate  that  sediment  attenuation  rate  is  linear  in  frequency 
requires  more  intensive  investigation,  especially  at  the  lower  frequencies. 

3.  The  M-S  model  performs  well  at  frequencies  above  1  kHz.  Below  500  Hz, 
where  bottom  and  subbottom  properties  are  important,  there  are  conflicting  reports  of  its 
performance.  At  long  ranges  and  low  frequencies,  for  instance,  the  Colossus  data  base  for 
the  Atlantic  Coast  shows  a  higher  loss  than  that  measured  in  South  China  Sea  areas. 

4.  Rogers’  computational  method  of  fitting  normal-mode  solutions  to  two  (or  three) 
propagation  loss  coefficients  has  some  advantages.  It  may  be  applied  to  sensitivity  anal¬ 
yses  where  parameters  are  changed  or  added.  It  also  furnishes  a  convenient  way  to  sum¬ 
marize  the  effects  of  treating  new  phenomena.  For  example,  the  Kuperman  and  Ingenito 
normal-mode  solution  to  the  boundary  scattering  problem  can  be  conveniently  compared 
with  reference  models. 

5.  Rogers’  model  and  the  extended  M-S  model  yield  close  values  of  propagation 
loss  when  the  same  inputs  are  used. 
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surface;  and  (3)  a  simplified  Rayleigh  two-fluid  model  of  the  bottom  for  sand 
or  mud  sediments.  With  a  few  free  parameters,  including  water  depth,  about 
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The  M-S  model  yields  good  predictions  when  applied  properly;  it  is  not  to 
be  used  for  all  environments. 

The  M-S  model  is  compared  with  Rogers*  semiempirical  model  based  on  numerical 
calculations  of  the  normal-mode  solutions  to  the  wave  equation  for  propaga¬ 
tion  in  shallow  water.  (Rogers'  model  for  transmission  loss  involves  a 
simple  range-dependent  fit  with  three  coefficients,  and  allows  a  convenient 
sensitivity  analysis  for  the  addition  or  changing  of  parameters.)  Differ¬ 
ences  in  predicted  propagation  losses  are  discussed. 

The  M-S  model  is  also  extended/ for  the  purpose  of  treating  arbitrary  nega¬ 
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